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We study the chiral phase transition at T ^ in the dual Ginzburg-Landau theory 
using the Schwinger-Dyson (SD) equation. In order to solve the SD equation at 
T ^ 0, we provide a new ansatz for the quark self-energy in the imaginary-time 
formalism. The recovery of the chiral symmetry is found at T c ~ lOOMeV with 
realistic parameters, which are set by reproducing the values of the string tension 
yf5 ~ 0.44GeV and the chiral condensate (qq) ~ -(250MeV) 3 at T = 0. 



The dual Ginzburg-Landau (DGL) theory is considered as an effective 
theory of QCD based on the dual Higgs mechanism in the abelian gauge 111. In 
the QCD-monopole condensed vacuum, its Lagrangian is described as 



4 V + 2 (n-d) 2 +n 2 m 



£dgl = —rt^+-f-A„ i ,. 2 , V 2 - A„+q(ip-e£H-m q )q, (1) 



where f M „ = d M A„ — d„A M and rt M is the direction of the Dirac string. (The 
notations are the same as those in Ref. 2.) In the Landau gauge, the gluon 
propagator in the QCD-monopole condensed vacuum {ms ^ 0) is given by 
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with X^y — -^2 s x ^ape \ ul sn a n' 1 k 13 k s . It is noted that the double pole factor 
( ra .fc)is in the second term, which is generated as a result of QCD-monopole 
condensation, produces the long-range coxrelation between quark and anti- 
quark due to its strong infrared singularity!!! The linear confinement potential 
is derived from this double pole within the quenched approximation in the 
static qq system El. However, this singularity causes the infrared divergence in 
the Schwinger-Dyson (SD) equation within thc-rainbow approximation, which 
corresponds to the quenched approximation Ela. In previous publications, we 
conjectured that the infrared screening effect for the double pole should be 
generated by the dynamics of light quarksBEl. 

Thus, we discuss the full gluon propagator 2? M „ in terms of the vacuum po- 
larization tensor n Ml/ as T)^ — — _D A1 A[n A<T ]2? cr „ . Owing to two conditions, 
Xftvk 1 ' = and X^n 11 = 0, we can sum over all the polarization diagrams and 
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obtain the full gluon propagator as 
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with n^fc) = [k^ky — k 2 g^ u )U s (k 2 ) and = m 2 B /{\ + n s ). We find that 
the infrared divergence at the quenched level is removed due to the vacuum 
polarization of dynamical quarks. Hereafter, we introduce the infrared cutoff a 
to the double pole in Eq.(||) as Em. We also take the angular average 

on the direction of the Dirac string&El for the gluon propagator Df£, with the 
infrared cutoff a as 
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The SD equation in the rainbow approximation and in the chiral limit is 
given by 

S^(P) =t>~i J -0^Q 2 rS q (k)^D;i( P -k) . (6) 

For the quark propagator, we can take S~ 1 (p) — z^ 1 (p 2 )[p' — M(p 2 )}; z(p 2 ) 
is the wave function renormalization and M(p 2 ) is the quark self-energy. The 
SD equation (^) is then decomposed into a pair of integral equations by taking 
the trace of S~ 1 (p) and ^S~ 1 (p). The coupled equations are obtained as 
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with fcp = p^ — k^ in the Euclidean metric after the Wick rotation in the 
fco-plane. Here, the Euclidean variables are simply denoted as p or k. In the 
case of ms = 0, z{p 2 ) is equal to one for all momenta in the Landau gauge, 
since the second term in the r.h.s. of Eq.(H) is equal to zero by taking the 
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angle integration □. We took therefore an passumption as z(p 2 ) = 1 for the 
sake of simplicity in previous publications Em. In addition, this assumption 
corresponds to the replacement of d{k 2 ) into d(max{p 2 , fc 2 }) in Eq.dq)tl 

In this paper, we solve numerically coupled integral equations (0) and (||) 
without above assumptfoiL As for the gauge coupling, we adopt the hybrid type 
of the running couplingEm, which is reduced to the perturbative-QCD running 
couplinga. We show in Fig.l the results on the wave function renormalization 
z(p 2 ) and the quark self-energy M(p 2 ) as functions of p 2 . This figure indicates 
that z(p 2 ) is almost equal to one in the QCD-monopole condensed vacuum 
(tub i= 0). Therefore, the assumption as z(p 2 ) = 1 seems reasonable. 
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Fig.l The wave function renormalization z(p 2 ) and the quark self-energy M(p 2 ) as 
functions of p 2 with e = 5.5, tub = 0.5GeV and a = 85MeV. The result with an 
approximation z(p 2 ) — 1 is shown by dashed curve. The unit is taken as Aqcd — 
200MeV 

Now, we formulate the finite-temperature SD equation using the imagi- 
nary time formalism to Eq.(Q) with z(p 2 ) — 1. Then, the quark self-energy 
M T (bj n ,p) depends not only on the three dimensional momentum p, but also 
on the Matsubara frequencies oj n - As a consequence, the resulting equation is 
very hard to solve even numerically. We propose the covariant-like ansatzB, 
in which a replacement is made for the quark self-energy at T ^ 0, instead of 
M(p 2 ) -> A/ T K,p), as 

M(p 2 ) -» M T (p 2 + col) = M T (p 2 ) (9) 

with p 2 = p 2 + w 2 and ui n = (2n + 1)ttT. It is noted that this ansatz guar- 
antees that the finite-temperature SD equation in the limit T — > is exactly 
reduced to the SD equation at T = 0. This fact is also confirmed by numerical 
calculations as shown in Fig. 2. The final form of the SD equation at T ^ is 
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derived as 
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where k 2 m = k 2 + p 2 — 2uy (p 2 — uj 2 )(k 2 — uj 2 n ) — 2uj m ui n - We then solve the 
SD equation by setting uo n = in the r.h.s. of Eq. (|Io[) I. We show in Fig. 2 
the quark self-energy M T (p 2 ) as a function of p 2 at finite temperature. No 
nontrivial solution is found in the high temperature regioni, T ~ HOMeV. In 
other words, the chiral symmetry is restored at high temperature. 

In summary, we have confirmed that the infrared singularity problem is 
avoided due to the screening effect of the quark vacuum polarization. By solv- 
ing numerically the coupled SD equations at T = 0, we have found that the 
wave function renormalization is almost equal to one in all momenta and the 
quark self-energy has non-zero value. We have solved the the SD equation 
at T ^ numerically with the covariant-like ansatz and found that the chi- 
ral symmetry is restored at high temperature; T c ~ lOOMeV with realistic 
parameters 
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Fig.2 The quark self-energy M T (p 2 ) as the function of p 2 at T = 0, 60 and lOOMeV. 
The same parameters are used as in Fig.l. 
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